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Abstract 

We generalize the Np = 2 Schwinger model on the lattice by adding a charged scalar 
field. In this so-called xU<p2 model the scalar field shields the fermion charge, and a neutral 
fermion, acquiring mass dynamically, is present in the spectrum. We study numerically 
the mass of this fermion at various large fixed values of the gauge coupling by varying the 
effective four-fermion coupling, and find an indication that its scaling behavior is the same 
as that of the fermion mass in the chiral Gross-Neveu model. This suggests that the xU(t>2 
model is in the same universality class as the Gross-Neveu model, and thus renormalizable 
and asymptotic free at arbitrary strong gauge coupling. 
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1 Introduction 



Strongly coupled gauge theories tend to break dynamically chiral symmetry, but fermions 
which acquire mass through this mechanism are usually confined, as it is the case in the 
Schwinger model or in the QCD. From the point of view of the electroweak symmetry 
breaking in, or beyond the standard model, a dynamical mass generation without the fermion 
confinement is of interest. Such a situation arises in a class of chiral symmetric strongly 
coupled gauge theories on the lattice, in which the gauge charge of the fermion acquiring 
mass dynamically is shielded by a scalar field of the same charge pp. The question is whether 
such models are renormalizable at strong gauge coupling, so that the lattice cutoff can be 
removed, and the resulting field theory might be applicable in continuum. 

In this work we consider such a lattice model with U(l) chiral symmetry and vectorlike 
U(l) gauge symmetry, the xU4'd model, in = 2 dimensions. It consists of a staggered 
fermion field a gauge field U G U(l) living on the lattice links of length a, and a complex 
scalar field with frozen length |0| = 1. The unconfined fermion field is F = (p'^x- There 
is no Yukawa coupling, and the mass amp of the fermion F arises dynamically. Because of 
the fermion doubling the fermion number is Np = 2. The xU<p2 model can be seen either as 
a generalization of the Schwinger model with Np = 2 by adding a charged scalar field, or as 
the 2D scalar QED with added fermions. Similar models have been studied nearly 20 years 
ago in the context of the instanton investigations [2]. 

The four dimensional xU 04 model, considered as a possible theory with dynamical mass 
generation of unconfined fermions has been investigated recently. In spite of some en- 
couraging results jHl 15 a- clarification of its renormalizability properties remains a difficult 
task. Here we demonstrate that in a simpler case, in 2 dimensions, these properties can 
be investigated with remarkable clarity. Within the limits of numerical accuracy, we find 
that the xU(l)2 model is renormalizable at strong gauge coupling g, because it belongs to the 
universality class of the two-dimensional chiral Gross-Neveu (GN2) model with Np = 2. 

The GN2 model is known to generate fermion mass dynamically at arbitrarily weak four- 
fermion coupling G, to be nonperturbatively renormalizable, and asymptotically free. One 
of its lattice regularizations is identical to the xU<p2 model in the limit of infinite gauge 
coupling, (3 = l/a^g"^ = 0. In this limit the four- fermion coupling G is an invertible function 
G{k) of the hopping parameter k, of the scalar field (or, equivalently, of its bare mass) 
in the xU<p2 model, with G{n) — s> as k ^ 00. The scaling behavior associated with 
the asymptotic freedom, and the continuum limit of the GN2 model, are thus obtained as 

These are extremely useful facts when the xU4>2 model is investigated at a finite gauge 
coupling, i.e., at nonvanishing (3. The idea is to compare the scaling behavior of the xU(t>2 
model with that of the GN2 model, as k grows at fixed /3 > 0. For this purpose we intro- 
duce an effective four-fermion coupling G{(3,G). This is now a coupling of the "composite" 
fermions F = (p^Xy ^ind thus characterizes the van der Waals forces arising from the fun- 
damental interactions between the fields. At /? = it coincides with the GN2 coupling, 
G{0,G) = G. For /? > 0, G is smaller than G, but depends on n similarly as G{k). For 
the comparison it is therefore convenient to use G instead of n as one argument of G. We 
investigate the scaling behavior of amp with decreasing G{/3, G) at various fixed (3 > 0, and 
compare it with the (3 = case. 

A determination of G{k, f3) directly by means of the four-point function of the composite 
field would be very expensive. Instead, we introduce it in an indirect way: we make the 
assumption that the scaling behavior of amp in the xU(t>2 model at /5 > is described by 
the truncated Schwinger- Dyson (SD) equations of the same structure as at /5 = in the GN2 
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model. The only change is the replacement of G by G{G,(3). Our numerical study of the 
xU4>2 model at /? > is mainly concerned with the verification of this assumption. 

The use of the SD equations serves further purposes. One is to provide an analytic 
framework in which the numerical results obtained at small but nonzero bare fermion mass 
amo can be related to the chiral limit case, tuq = 0, we are actually interested in. Simulations 
at tuq = are too time consuming on larger lattices, and the scaling behavior has to be 
studied at arriQ > by varying both G and amg. The SD equations suggest how to vary 
both parameters simultaneously, approaching the critical point, G = aniQ = 0, and how to 
extrapolate to the chiral limit mo = 0. We test this strategy carefully at /5 = 0, exploiting 
the knowledge of the chiral limit there, and find it remarkably successful. Small deviations, 
not accounted for by the SD equations within the used truncation scheme, and present at 
very small arriQ, do not influence the correct scaling behavior following from the asymptotic 
freedom. 

Further benefit of the use of SD equations is the control of the finite size effects. Following 
Ref. we solve these equations on lattices of the same sizes and boundary conditions as 
those on which numerical data are obtained. This allows a suitable choice of the data points 
in the parameter space, as well as an extrapolation to the infinite volume limit. 

Finally, the SD equations allow us to describe the preasymptotic behavior of the data, 
obtained for correlation lengths limited by the lattice size, and to infer the genuine scaling 
behavior for diverging correlation length. At j3 = 0, we have verified that when the SD 
equations are used for the extrapolation, then from those data the correct scaling behavior 
with G — > is obtained. 

We find that in the xU<p2 model with (3 < 1 the preasymptotic behavior, and thus 
presumably also the scaling behavior of the fermion mass amp is described quite well by the 
same SD equations as in the chiral GN2 model, now with G{(3, G) replacing G. The same 
scaling behavior as in the GN2 model is indicated, with G ^ as /t ^ 00. 

On the basis of this evidence we suggest that the xU 02 model belongs, along the critical 
line at K = cx), to the same universality class as the GN2 model, and is thus nonperturbatively 
renormalizable. The xU4'2 model is therefore an example of a renormalizable quantum field 
theory in which the fermion mass is generated dynamically at strong gauge coupling by the 
shielded gauge mechanism suggested in Ref. p. 

In the following section we define the xU(f)2 model and explain its relationship to the GN2 
model at /5 = 0. In Sec. 3 the phase diagram is described. The SD equations, and also the 
determination of the effective four-fermion coupling G(/3, G) by means of their inversion, are 
discussed in Sec. 4. In the Sec. 5 we test these equations, and the accuracy of determination 
of G at /9 = 0. In Sec. 6 we then apply the same method of analysis to the data for amp 
at < /3 < 1, demonstrate the applicability of the same SD equations as at /? = 0, and 
determine G(/3, G). The scaling behavior of the data is illustrated in Sec. 7. We conclude in 
Sec. 8 with some remarks about the meaning of our results. 

2 The xU(t>2 model and its GN2 limit 

The action of the xU(t>2 model consists of three parts: 

SxU(j) = Sy. + Su + S^, (2-1) 

where 

1 ^ 

•^X ^ y Xx ^ y Vx^{Ux^^Xx+n Ux_^^^Xx—^) y ^ XxXx ) 

^ X fl=l X 
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Figure 1: Phase diagram of the xU(f>2 model for ttiq = 0. The weh understood limit 
cases are described. The dashed line indicates a possible topological phase transition. The 
fermion mass amp is nonzero everywhere except on the bold-marked boundaries. 



Su = f3j2i^-ReUp) 



X 11=1 

Here, x and x + fi denote lattice sites and their nearest neighbors, respectively, and x is a 
staggered fermion field, r/^.^ being the usual phase factors. The link variables U^^^ G U(l) 
represent a compact abehan gauge field with coupling P = and is a complex scalar 

field with the constraint l^^^l = 1. The fermion and scalar fields have the same unit charge. 
The bare fermion mass mg is introduced for technical reasons, and we are interested in the 
limit rriQ = 0, where the action has a global U(l) chiral symmetry. The phase diagram of 
this model is shown in Fig. ^ 

In the limit case (3 = 0, one can perform the Lee-Shrock transformation [Hj, in which the 
scalar and gauge fields are integrated out and a four-fermion term appears. This leads to 
the action 



4f 



EE G XxXxXx+fiXx+fj, 

X ^=1 \ 
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:Vxfi 



with 



and 



G = G{k) 



XxXx+fj. Xx+fiXx 



amp 

+ 2^ XxXx 

X 



4r2 



/i(2/t) 



(2.2) 



(2.3) 



(2.4) 



For mo = the transformed action ()2.2j) is that of the chiral GN2 model in a certain 
lattice regularization. Within the uncertainty of interpretation of a continuum limit of 
staggered fermions with a strong gauge coupling it can be possibly interpreted also as a 
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G(/c) = (l-r'^)/4r^ r(/c) =Ii(2/c)/lo(2/f ) 




Figure 2: Relation between the four fermion coupling G and k at /? = 0. 



lattice formulation of the Thirring model The GN2 model has a critical point at G = 0, 
where the fermion mass vanishes like 



amp oc e 



TT 



(2.5) 



The four fermion coupling G{k) is a function of k shown in Fig. |21 From the point 
of view of the scaling behavior ()2.5|) . the use of G{k) as a parameter instead of k is very 
convenient, and we therefore adopt such a reparametrization even at /? > 0. There G is not 
a four-fermion coupling any more, it only replaces the hopping parameter k according to 
eq. ()2.3|) . Also r, eq. ()2.4|) . will be understood as a function of G from now on, satisfying 
r(G) ^ 1 as G 0. 

We note that the bare mass amg/r in the action ()2.2j) of the GN2 model is different from 
that of the xU(t'2 model ()2.1|1 . as the field x has been rescaled by ^/r in the course of the 
Lee-Shrock transformation. We use amo of the xU<p2 model, and the bare mass in the SD 
equations for the GN2 model is therefore slightly G-dependent at fixed amo. 

We perform the hybrid Monte Carlo simulations on the V = L'^ lattices, with periodic and 
antiperiodic boundary conditions in first and second (/i = t) directions, respectively. The 
fermion mass amp is obtained from the gauge invariant propagator {4>l;Xx<PyXy) by inverting 
the fermion matrix M and calculating 



x,y 



(2.6) 



k=0 



Here tta denotes the usual momentum shifts in the Brillouin zone. The mass amp, and also 
the fermion renormalization constant Z result from the fit in momentum space to 

—4i sin kf 



TYTtP{h) = Z 



sin^ kt + (arriF) 



(2.7) 



where Fj is the Golterman-Smit matrix {T^=t)AB- This procedure is chiral invariant. 
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3 Other limit cases and the phase diagram 



The schematic phase diagram for mo = is shown in Fig. ^ At k = 0, when the scalar 
field is decoupled, the xU(t)2 model with mo = reduces to the lattice regularized Schwinger 
model id = 2 QED) with Np = 2. As is well known, this model is anomaly free, confines 
fermions at all /3 < oo, and possesses only massive bosonic states. Thus amp — > cxd as k — 0. 
The continuum limit is expected at the UV fixed point at /? = cxd, and there is no critical 
point at /? < oo. 

At K = oo, both the scalar and the gauge fields are frozen0 As is seen in the unitary 
gauge, the fermion field is decoupled, and amp = am^. At amo = 0, the line k = oo is thus 
a critical line at which the fermion mass amp vanishes. The /? = point of this line is the 
critical point of the GN2 model. 

In the limit (3 = 00, when the gauge coupling vanishes, the fermion field is decoupled 
again, and amp = amo. Thus, in the chiral limit, the line /? = 00 is a critical line with 
vanishing fermion mass, too. The scalar field variables can be seen as spin variables, and the 
corresponding two-dimensional XY2 model is known to have a topological phase transition 
at K = ~ 0.56. 

Let us now consider the inside of the phase diagram in Fig. H Old investigations suggest 
that the model possesses a massive fermion in some parameter region accessible to the dilute 
instanton gas ^2]. As follows from the convergence of the strong coupling expansion, at 
small nonvanishing (3, the model should have the same properties as at /? = 0. This implies 
analyticity and nonvanishing amp for k < 00 at small /3. 

The fate of the topological transition at /5 = 00, as /5 gets finite, is not completely clear. 
We have investigated numerically the spectrum of the model in the vicinity of the dashed 
line shown in Fig. ^ This line is observable at large /3 as a shallow dip in the masses of the 
scalar and vector bosons, which can be constructed from the gauge invariant products of the 
type (pl.Ux,^4>x+iJ,- But amp shows no sensitivity when the dashed line is crossed, and we have 
found no state, neither bosonic nor fermionic, indicating a vanishing of the mass in lattice 
units on this line at /? < 00. The dynamical fermion mass generation at /5 < 00 is thus 
not infiuenced by the remnant of the Kosterlitz-Thoules transition. Presumably, a critical 
behavior on this remnant, if any, appears only in some topologically nontrivial observables 
jH]. The mass amp stays finite, and the fermion F gets infinitely heavy in physical units in 
any conceivable continuum limit taken on this line at /3 < 00. 

We have checked that there is no indication of any other phase transition, not even of 
some change of behavior of some local observable, anywhere else in the phase diagram of the 
xU4>2 model. Our data thus indicate that the fermion mass amp is nonzero for any finite 
f3 and k. It decreases when any of these parameters increases and one of the boundaries 
bold-marked in Fig. ^ is approached. As in the GN2 model, the fermion mass generation 
takes place without the spontaneous chiral symmetry breaking, which is forbidden by the 
Mermin- Wagner- Coleman Theorem. A continuum limit taken on the line (3 = 00 might lead 
to an interesting generalization of the Schwinger model in the continuum. In this work, we 
concentrate on the scaling behavior of amp when the line k = 00 is approached. 



^ This is meant in a sense avoiding the Mermin- Wagner-Coleman theorem: the Hmit k ^ 00 is taken 
before an external "magnetic" field, required for a definition of the spontaneous symmetry breaking, is 
switched off. We thank G. Roepstorff for a discussion on this point. 
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4 Schwinger-Dyson equations for the GN2 model 

The SD equations for the fermion propagator in a four-fermion theory, truncated after 0{G), 
can be represented graphically as 



(±) 



+ 



On a finite lattice they read (see e. g. 0): 



(4.1) 



N=^ + ^T N 

r{G) V ^Y.uF^{^HKa)f + N^' 

F, = 1 + ^E ^'"^^^^y • (4.3) 

These are three coupled equations for and which we solve numerically. Then 

amp = N/Ft , (4.4) 

Z=^ (4.5) 

rFt ^ ' 

are determined. If the term (b) in ()4.1|1 is neglected, one obtains the gap equation with 
F^ = 1 and amp = N. 

In the infinite volume, when Fi = F2 = F, and for small G, the approximate analytic 
solution of the SD equations is 

amo 8GN ^ /2N\ ,^ , 

^ = ;;Tln^ , 4.6 



F = - + -VTTG. (4.7) 

For mo = 0, one obtains the scaling behavior ()2.5j) as G — 0. 

The idea of the combined limit G — and arriQ — * is to make aruQ/r a function of G 
in such a way that eq. ()4.6|) is solvable, and that amo(G) — as G ^ 0. We choose 

= (I - sf-iy^T , (4.8) 

using the approximate solution ()4.7p for F. Here s is a free parameter obeying < s < 1. 
For this choice of amo the solution of eq. ()4.6|1 is 

N = —e~^^, (4.9) 
2 
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and thus 

N n 



amp = — = -e 8G . (4.10) 

r 2 



The asymptotic scahng law along the lines of constant s is 



TVS 

amp = —e 8G , (4-11) 
amp 1 



arriQ 1 — s ' 



(4.12) 



where we have used rF ^ 1 as G ^ 0. Obviously, the ratio mo/rriF does not vanish except 
if s = 1. For s < 1, a continuum GN2 model with nonvanishing bare mass is obtained. We 
therefore distinguish between the chiral limit, niQ = for any a, and the statement arrio = 0, 
which at a critical point with amp = allows a = 0, mo/mp 7^ 0. 

It may appear strange that the scaling behavior 1)4.111) depends on s, i. e. on the path, 
though the same critical point is approached. But this is similar to the scaling behavior in 
a magnetic system described by the equation of state 

h = M'f (-^) , (4.13) 



when the external field h (corresponding to amo/r) and the reduced temperature t (corre- 
sponding to G) are varied simultaneously. For example, if 

h = c-tP, P<I35, (4.14) 

the scaling law is p-dependent, 

M oc fp'^ . (4.15) 

In our case the scaling is described by an essential singularity instead of the power law ()4.15p . 
Needless to say, the SD equations reproduce correctly the first coefficient of the /9-function 
of the GN2 model with nonvanishing bare mass since the one loop contribution is taken into 
account correctly. 

The truncated SD equations cannot be expected to describe the data correctly at very 
small amo, because the contributions of the neglected fermion loops increase with decreasing 
arriQ. Then only the exponential scaling behavior, eq. ()4.1H) . but not the value of the constant 
prefactor in ()4.1H) . and thus of the ratio ()4.12|) . are predicted correctly. We shall take this 
discrepancy into account by relying in our determination of G in Sec. IHl on the results at 
larger values of amo (mostly amo = 0.4). 

An important step in our data analysis at /3 > is the use of the inverse SD equations. 
For a chosen /5, s and amo, the fermion mass amp, obtained on a certain lattice, is inserted 
into the corresponding SD equations ()4.2|) . ()4.3p by using eq. ()4.4p . The four fermion coupling 
is considered as a free parameter F. It is determined by solving the three equations ()4.2j) . 
()4.3j) . taken as equations determining F and F^: 

^ p ^ Q"^o I 1£ ampFt 

"""^^ ' ~ r{V) F V Fl (siniKa))' + {ampF.f ' ^ ' ^ 

F, = 1 + ^V F^jsHKa)? _ 

y k EuF^ (sin(M) + {ampF.y ^ 
The resulting values of F can, in principle, depend on all the parameters, F = F(/3, G, amo, V)- 



8 



5 Test of the SD equations in the GN2 limit 



The aim of this section is to investigate how the known scahng behavior ()2.5j) of the GN2 
model can be confirmed by the data for amp obtained in numerical simulation of the xU(t>2 
model at /5 = 0. The fermion matrix inversion has turned out to be very slow at rriQ = 
on large lattices (> 32^), and the simulations have to be performed at finite arrio. We have 
made several attempts to extrapolate to ttiq = the values of amp obtained for several 
amo at fixed G. However, using e. g. a power law we failed to reproduce reliably the values 
obtained by long simulations directly at mo = 0. 

Therefore we have adopted the strategy of the combined approach to the critical point 
G = arriQ = 0, in which G and arriQ vary simultaneously, choosing the paths s = const, 
eq. ()4.8p . suggested by the SD equations. The chosen values of s are s = 0.2, 0.3, 0.7. The 
simulations have been performed at the values of amo and G satisfying the relation ()4.8j] . 
and chosen such that the value of amp predicted by the SD equations satisfies 

1 L , , 

< T • (5-1) 



amp 4 

This restriction turned out to be necessary in order to avoid a finite size dependence both 
of the predicted and measured values of amp. 

In Fig.Elwe show the semi-logarithmic plot of the data for amp against l/G, and compare 
them with the prediction of the SD equations ()4.2|) and ()4.3|1 on the 32^ and 64^ lattices. 
It demonstrates that the SD equations describe the data for smaller s very well, and for 
larger s, s = 0.6, 0.7, still quite well, though for these s the bare mass is very small. But one 
cannot see that the data is still far from the asymptotic scaling behavior ()4.1H) . and that 
the seemingly linear decrease does not have the right slope. Also the size of the error bars 
is barely visible on the logarithmic scale. 

Therefore, in Fig.|3]we show the data as ratios amp /amo plotted on the linear scale. The 
curves are the SD predictions for this ratio on the 32^ and 64^ lattices (full lines) and 1024^ 
lattice (dotted lines), whereas the dashed horizontal lines represent the expected asymptotic 
scaling behavior, eq. ()4.12j) . The benefits of the SD equations become manifest: For smaller 
s the agreement with the data is within the tiny but now visible error bars. With increasing 
s the agreement gets worse, as expected for decreasing amo. Nevertheless, the SD equations 
still reproduce the l/G- dependence qualitatively. They predict the onset of finite size effects, 
manifested by a downward bend of the curves. These effects set on at larger l/G on larger 
lattices. We have checked that the data not satisfying the restriction ()5.H) behave in a 
similar way. Furthermore, one can observe the now apparent difference from the asymptotic 
scaling behavior ()4.12j) . with an indication how this behavior would be slowly achieved on 
huge lattices (dotted lines). We find that the SD equations describe the data for s < 0.5 
long before the asymptotic scaling sets on. This allows us to extrapolate the data obtained 
on achievable lattices to large l/G, and to interpret them as an evidence for the asymptotic 
scaling ()4.11|) - ()4.12|) . expected for the GN2 model with nonzero bare mass. 

We expect that the discrepancies observed at s = 0.6 and 0.7 are due to the truncation, 
and thus do not indicate any deviation from the asymptotic scaling. This is further supported 
by the observation that the agreement between the data and the SD equations is improved 
even at larger s, if we plot the data for amp/am^Z , with Z measured by means of the 
relation ()2.7p . and compare them with the SD equation using ()4.5p (see Fig. Because the 
measured values of Z are consistent with Z ^ 1 for G ^ 0, as follows also from eq. ()4.5|1 . the 
asymptotic behavior is not changed, and we now see how it is approached even for larger s. 

As a preparation for the studies at /? > 0, it is illustrative to test, at /? = 0, also the 
inversion of the SD equations according to eqs. ()4.17|) and ()4.16|) . The above data for amp 
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Figure 3: The data for amp at /? = for various fixed s compared with the predictions of 
the SD equations. 
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Figure 4: The ratio anip/aniQ at /? = for various fixed s compared with the SD equations 
on 32^ (s < 0.5) and 64^ {s > 0.5) lattices (full hues), and 1024^ lattice (dotted lines). The 
dashed horizontal lines represent the expected asymptotic behavior. 
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13 = 0.0 




has been used as input to these equations, and T{0,G,amo,V) has been obtained by their 
numerical solution. The results are shown in Fig. IHl where they are compared with the 
expected value of the coupling G{0, G) = G. We observe that for a rather large bare mass, 
arrio = 0.4, the agreement is excellent. The fit r(0, G, 0.4, V) = cG^ gives, for V = 32^ and 
64^, results consistent with c = p = 1 (see table 1). 

For smaller amo, deviations of F from the true coupling G are observable. Of course, 
these deviations are of the same origin as those in Figs. El and |3 they reflect the inaccuracy of 
the SD equations for larger s, and thus for smaller arriQ. For a quantitative comparison with 
the SD equations it is therefore advantageous to perform simulations at not too small amg. 
Nevertheless, even for smaller amo the deviations remain small. The degree of agreement 
between F and G at f5 = can serve at > as an estimate for the accuracy, with which one 
can determine the effective coupling G(/5, G) from the values of F(/5, G, amo, V) obtained by 
the inversion of the SD equations. 

6 Effective four-fermion coupling at > 

The experience gained at /5 = makes it clear that also at /5 > the simulations have to be 
performed at amo > 0, and that on lattices of affordable sizes there is no chance to observe 
the asymptotic scaling directly. Therefore, we investigate the scaling behavior of the fermion 
mass amp at /5 > by means of the following strategy: 

1. The conjecture that the xU<p2 model at /5 > belongs to the same universality class 
as the GN2 model suggests the use of the same SD equations as at /3 = for the 
description of the nonasymptotic data for amp- The only foreseen difference is the 
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Figure 6: The values of r(/5, G, amo, ^) obtained by the inversion of the SD equations at 
/? = for V = 64^. The fuh Une representing the fit by means of eq. (|6.2|) to the arriQ = 0.4 
data nearly coincides with the diagonal (dashed) line. 



value of the effective four fermion coupling G{l3,G), which is not any more a known 
function ()2.3|) of k, but has to be determined from the data. A similar idea was very 
successful in the earlier study of the universality class of the two-dimensional Yukawa 
model p. 

2. For each value of amp, obtained at some (3, G and amo on a lattice of volume V, we 
invert the SD equations ()4.1f)j) and ()4.17|1 . obtaining T{P,G,amQ,V). 

3. Of course, in principle the effective coupling should be independent of amo and V. 
Because of the limited accuracy of the truncated SD equations, some dependence of 
r{l3,G,amo,V) on amo remains, however. On the basis of the observation that at 
P = the values of T and G are consistent for amo = 0.4, we assume that also at /? > 
the values of G) can be obtained for this amo, and define 

G{(3,G) =T{f3,G,amo = OA,V). (6.1) 

This effective four-fermion coupling is V^-independent, provided V is sufficiently large. 

4. At each P, we determine the G-dependence of G{P, G), which turns out to be consistent 
with a power law 

G{p,G) = c{j3)GP^^\ (6.2) 

5. If G{(3, G) were known prior to the numerical calculations of amp at /5 > 0, one would 
choose the data points so that in the {amo, G) plane they lie on lines similar to ()4.8p . 
with G replaced by G{[3,G). The same simple comparison with the expected scaling 
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Figure 7: The values of T{f3, G, amo, V) at (3 = 0.3 and /? = 1.0. The full lines are fits by 
means of eq. (|6.2j) to the arrio = 0.4 data. They represent the obtained effective coupling 
G{P,G). 



behavior ()4.11|) . now with G{l3,G), would then be possible. Without this knowledge, 
we have decided to acquire the data on the same lines of constant s, eq. ()4.8|) . as at 
/3 = 0. As will be explained in the next section, one can then recalculate the predictions 
of the SD equations with G{P, G) to the lines s = const., and observe the approach to 
the scaling behavior, though these lines are now less suitable for this purpose. 

We have determied amp at /3 = 0.2, 0.3, 0.5, 0.7, 1.0 on the 32^ lattice for s = 0.2, 0.3, 0.4, 
0.5, and at /3 = 0.3 also on the 64^ lattice for s = 0.4, 0.5, 0.6, 0.7. The intervals of the G and 
amo values have been chosen such that the fermion mass is consistent with the requirement 
()5.H) . As expected, the fermion mass amp decreases with increasing /3 at constant G, because 
of the decreasing gauge coupling. 

For (3 > 1 the mass amp is measurable only in a small G-interval below or in the vicinity 
of the dashed line in Fig. ^ Nevertheless, one can see that amp is insensitive to this line. 

In Fig. 13 we show the results for r(/3, G, arriQ, V) at P = 0.3 and (3 = 1. They cluster for 
each (5 along, or slightly below a simple curve. This curve is the power law fit ()fi.2j) to the 
values of G(/?,G), determined according to eq. ()fi.l|l at arriQ = 0.4. The results for smaller 
arriQ do not deviate from this curve more than in the (3 = case (Fig. IH)). For other (3 values 
the results are very similar. 

Fits by means of the power law ()6.2|1 describe the data at arriQ = 0.4 very well. In table [T] 
we present the results for the fit parameters c{(3) and p{(3). The indicated errors come from 
the Minuit fit for G{(3, G) by means of eq. ()().2|1 . As the values of G{(3, G) result from a data 
analysis in several steps, these errors are too naive. A more realistic error estimate might 
be the difference between the results on 32^ and 64^ lattices at (3 = 0, where we known 
that these results should be consistent. This suggests the errors of sizes Ac/c ~ 0.05 and 
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322 


642 


c 


P 


c 


P 


0.0 


1.038(7) 


1.016(4) 


0.99(4) 


0.99(2) 


0.2 


0.767(9) 


0.940(5) 






0.3 


0.675(5) 


0.913(3) 


0.62(1) 


0.88(1) 


0.5 


0.512(6) 


0.844(4) 






0.7 


0.422(3) 


0.806(2) 






1.0 


0.321(2) 


0.747(2) 







Table 1: The values of the parameters determining the effective coupling G(/3, G) by means 
of eq. (|6.2|) . See the text for a discussion of the errors. 



Ap/p ~ 0.03. 

These results demonstrate that the introduction of the effective four fermion coupling 
G{(3, G) by means of the SD equations is sensible. In other words, one can find such a 
function G{l3,G), that the data is consistent with the SD equations when G{[3.G) is used 
as a coupling. Assuming now that this is true also at smaller G and amo, i. e. beyond the 
intervals we could investigate, we can deduce the scaling behavior of amp with G(/3, G) from 
these SD equations. 

The power law dependence fl6.2j] of G'(/3, G) on G is questionable for very small G because 
of the singular behavior at G = 0, and probably should not be trusted there. K.-I. Kondo 
has investigated in the continuum a model quite similar to the xU4>2 model, both in d = A 
j7] and d = 2 (TU]. Recently, he has obtained in d = 2, by solving the SD equations for 
the full xU<p2 model, the effective four-fermion coupling G as a series in G ^U]- This would 
contradict ()6.2|1 . We have checked that our results for G can be described by a polynomial 
in G, too, but the coefficients of the higher powers of G are large and unstable, and such an 
analytic description is therefore of little use. 



7 Scaling behavior at /3 > 

According to the results presented in the previous section, the SD equations with G replaced 
by G{P,G) describe the data for amp for larger amo values. These equations predict also 
the behavior of amp /amo as G and amg approach zero along the lines 

1 ~ ^ ^ ttF ^ 

— amo(G, s) = (1 — s)——e 8G , s = const. , (7.1) 



r(G) ' ^ ^ '2 

in the (amo,G) plane. Along these lines the asymptotic scaling behavior is analogous to 



dnu, 

ITS 

TT amp 1 

amp = —e oLr , = z ■ (7.2) 

2 amQ 1 — s 

This is the same scaling behavior as in the GN2 model with amo > 0, when G is replaced by 
G. From eqs. ()2.3|) . ()2.4|1 and ()6.2|1 it follows that G{[3,G) — > as k — > 00. The successful 
analysis of the data for amp at /? > by means of the SD equations with G(/3, G) is thus an 
indication that the xU(t>2 model at /? > belongs to the same universality class as the GN2 
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model, when the critical line k = cxd is approached. In the chiral limit, the scaling behavior 
predicted by the SD equations with G{I3, G) is 



amp oc e 8G = e ScG^ , (7.3) 

to be compared with ()2.5p at /9 = 0. 

Of course, our data with rather large amo and amp do not show directly this asymptotic 
scaling, but only an approach to it, as predicted by the SD equations. To illustrate this, in 
Fig. IHlwe compare the data obtained on lines s = const, ()4.8|1 . with the predictions for these 
lines of the SD equations with G{(3, G). As in Fig. El we plot the ratio amp/amo, and the 
full lines are the numerical solutions of the SD equations with G{P, G). As at /5 = 0, Fig. 
we observe a very good agreement for smaller s (larger amg), whereas at larger s at least a 
qualitative behavior of the data is reproduced. 

The data at /5 > has been taken on the lines ()4.8j) . which are not as convenient as the 
lines (jT.lj) would be. Nevertheless, one can obtain the asymptotic behavior on the lines ()4.8|1 
as follows: comparing eqs. ()4.8p and ()7.1|) we eliminate amo, obtaining a relation between 
the pairs of parameters (G, s) and (G,/?), 

(l-s)e~8G^ = (l-5)e~8G^ (7.4) 
Using this relation in eq. ()7.2j) we obtain for small G (i.e. setting r = = 1) 

TT 

where s is now understood function s{f3, G, s). 

This function s(/3, G, s) can be determined at each f3 from the now known G-dependence 
()(i.2j) of G{P, G). The results at /5 = 0.3 for various fixed values of am^ are shown in Fig. IHl 
We observe that, for small s, the values of s are nearly independent of amo, and thus of 
G. This explains why the G-dependence of the data in Fig. |H1 looks very similar to the 
G-dependence of the data at /3 = 0, fig. up to a - nearly constant - factor (1 — s)/(l — s). 
This similarity is easily observable JI], as no determination of s(G, s) is required, and thus 
constitutes a simpler comparison of amp in the xU(t>2 model at /? > with the predictions 
of the SD equations with G. 

The observed agreement between the nonasymptotic data and the SD equations of the 
GN2 model with G suggests that also the asymptotic scaling behavior of the x^02 model at 
/3 > is described by these equations, and the model thus belongs to the same universality 
class. But, of course, this is only a conjecture, as the region investigated is limited by the 
applicability of the truncated SD equations, as well as by the constraints of a numerical 
approach. 



8 Conclusions and discussion 

By introducing an effective four-fermion coupling G(/3, G) in the xU<p2 model at /5 > we 
have found that - within the limits of a numerical approach - this model is described by the 
same SD equations as the chiral GN2 model with coupling G replaced by G. The scaling 
behavior, when G — >^ 0, is thus the same as in the GN2 model in the G — limit. This is 
an indication that the xU4>2 model belongs to the universality class of the GN2 model. 
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Figure 9: The values of s{f3,G,s) at (3 = 0.3 for fixed omo, determining G by using 
eq. (|4.8|) . The curves are independent of the lattice volume. The data are from Ref. lUj. 



From this we tentatively conclude that the xU4>2 model is nonperturbatively renormaliz- 
able and, though defined on the lattice, possesses a well defined continuum limit. It is thus 
an example of a strongly coupled gauge theory with continuous chiral symmetry, in which 
the fermion mass is generated dynamically and the massive fermions are not confined. At 
least in two dimensions the shielded gauge mechanism of fermion mass generation, suggested 
in Ref. P, exists. The xU<p2 model nicely illustrates this mechanism. 

On the other hand, the xU(t>2 model is presumably too simple to provide useful hints how 
to approach the d = A case. Our results suggest that the complexity of the model, 
and the composite structure of the fermion F = 0"^%, are effects not surviving the continuum 
limit. The possibility of describing the scaling properties in terms of the van der Waals force, 
represented by the effective four-fermion coupling, shows that in the renormalized theory the 
inner structure of the fermion, acquiring its mass dynamically, is irrelevant. 

The study of the renormalizability properties of the strongly coupled xU(f)2 model has 
been made relatively easy in = 2 by its neighborhood to the well understood GN2 model. 
Still, the applicability of the SD equations long before the onset of asymptotic scaling has 
been a surprise, even in the GN2 model. The use of these equations has been crucial, since 
an asymptotic scaling behavior is evidently not obtainable in numerical simulations, and a 
method of extrapolation to the scaling limit, and to the limit of chiral symmetry, is required. 

A plausible explanation why the xU 02 model is well described by the SD equations of the 
GN2 model might be as follows: Integrating out the gauge and scalar field, one ends up with 
a pure fermionic theory with many multi-fermion couplings. These fermions correspond to 
our fermions F, as in the (3 = limit. The universality of multi-fermion couplings in d = 2, 
observed e.g. in the studies of the d = 2 Yukawa model 0, suggests that the four-fermion 
coupling is dominant and sufficient to describe the data. 
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Provided that the xU 4>2 model at (3 > really belongs to the GN2 universality class, the 
most interesting question left is the dependence of the effective four-fermion coupling G on 
G and /3. We hope that our data for G{(3, G) will stimulate its theoretical investigation. 
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